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ABSTRACT 


The developing flow in the hydrodynamic entrance region of an equi- 
lateral triangular channel was investigated. Using a fully developed 
velocity profile from Knudsen and Katz (3) as a boundary condition, the 
equations were solved numerically employing the method of Chorin (1). 

The resulting velocity profiles and pressure drop were employed 
in calculating the local friction factor in the hydrodynamic entrance 


region. 
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I, INTRODUCTION 


The flow characteristics of non-circular channels have become 
of increasing importance in the analysis of compact heat exchangers, 
Nevertheless, heat transfer analyses of the triangular entrance region 
have been limited to approximations which postulate fully established 
velocity profiles (Graetz Approximation). Such solutions have been 
obtained by Lu and Miller (4) and Tao (8) for equilateral triangles, 

Using Sparrow's (7) results for the fully established velocity 
profiles in an isosceles triangular duct, McComas (5) has estimated 
values of entrance lengths and friction factors for such flows, 

In order to complete the analysis of heat transfer in the entrance 
region of triangular ducts, it is necessary to obtain a hydrodynamic 
solution for the developing velocity profiles in that region. It is 
therefore the purpose of the present work to provide a solution of 
the equations of motion in the hydrodynamic entrance region of an 


equilateral triangular duct, 
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II. ANALYSIS 


In an attempt to find a convergent solution to the equations of 
motion, several explicit methods were tried without success, Subse- 
quently a relaxation method employing an artificial density change, 
suggested by Chorin (1), was adopted. This method essentially alters 
the solution of the equations of motion from an initial value problem 
to a boundary value problem, necessitating the use of entrance length 
and the fully developed velocity profile as boundary conditions, 

A. Governing Equations 

The governing equations for the flow of an incompressible viscous 

fluid (Navier Stokes Equations) may be written: 


x —- Momentum: 


Su du du, yeu. i! <P 2 

oe an Ne 3 ayy, (1) 
y - Momentum: 

Wu wa a ur ov , . 

ot 3x 3y 34 § dy (2) 
Z - Momentum: 

re BM dur LAP) tur 

dt dX JY dy § 34 (3) 
Continuity: 

dU ge yw L 

yee ‘i 
ye” 4  s8 (4) 


12 





Introducing the dimensionless variables, 


us e 7. W 
ati ve Us W> ° 


Cl 


— eS 
Seige me 
P=? aa, makes 


equations 1-4 become: 


a yu a0) wee 2 
+R |o VN | 3% 


SY ay 


SW ab z. 
ew fe — St + VOW 
WaRlu gy a 32 


WwW. UW. 
x 3 ua ¥: 


Verluy Woy W, a Cer lte, 


B. Boundary Conditions 


A condition of uniform velocity is assumed at the entrance and 
no-slip conditions are imposed at the walls; thus: 
1. at the entrance:* 
MON 7) = 
V (0,Y,7,t) = 0 
WialO/ Nez, t) 7 0 


2. at the walls: 


Ut) = Va)=-W)= 0 


*The prime has been dropped and dimensionless variables are 
understood. 
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(5) 


(6) 


(7) 


(8) 


(9) 


CC. Antaidaeial Density 


Because the problem is one of incompressible flow, no density 


change can take place. One can introduce an auxiliary continuity 


equation, os f 35) ‘a OW | SV 


sts ee =O 
ay cee, Se oY 


(9a) 


in which changes of the pseudo density given by 


BES (10) 
can be used to test for convergence at each iteration process, 
Thus, we have introduced a set of auxiliary parameters, 5 and § |, 
which are analogous to an artificial density and an artificial com- 
pressibility. 

By introducing the artificial equation of state (equation 10) 

and using t’ as an auxiliary parameter analogous to time in compressible 
flow, the initial value problem can be altered to a boundary value 
problem. Such a change, however, requires the use of a fully estab- 
lished solution as an additional boundary condition. An exact solution 
for the steady, fully established profile has been obtained by Knudsen 


and Katz (3): 
UY,2) = -26R[¥(7+F)-g(7- 32F )- S| 


ye. oy a _ . 
where: 2=Y BEY and: p= 2-95 o The friction factor f for fully 


(11) 


developed flow is given by Kays (2) as 13.33. This number serves to 
determine the value of the fully established pressure gradient. 
D. Finite Difference Equations 

Equations 6=9 were put into finite difference form using central 
differences, except for the second derivatives which were differenced 


using the Dufort=-Frankel pattern, This latter expression is of the form: 


14 








‘qt Be Uy tx) +U'%-)- eee 
S55 (Sie 


where UN(x+) denotes U(N,X+4xX,Y,Z), UN+l is U(N+1,X,Y,Z) and N 


(12) 


denotes the Nth relaxation, The resulting equations were solved for 


the N+l terms and the resulting computing equations are of the form: 
= j1/( (ee oe bat )} x 1-5 or g | Ser - Fu- | 
- REE [Uke Tee - Dia deo] — BS [UdeVEn - Us VE] 
; Re [Uo We -Ua-1We | =r Zen, [Ue +Ue>- ie 
+ ae 2 | Ua) + UG - i a fn se |Ue4 Ue a “y] Cis) 


yun: (ee ee re 2a x os <[ 9" tye) = Fv) 
~ R= ax {v" ee Ve k-) Uk- | - Re * [Veer Vern - VeVor] 
= Re VenWes - Veo We] ae Ee, [Wee #V- ie | 

a ae Vern OV ie wee | + tet Nes rVeo-ve wey 


H L 
Ww as 1+ BR + A Oe )} we ae § [SED - SE | 
Re 7 [Wee Ube Wk Ue-» | _ [oe LY Y fw 4} Vw We Vero | 
2 [Way Weo We swe] + Far [wee Woo Ww" | 
(15) 


+ JF ye a Q- 
+e [Wen +Wey-y We J+ 2 zat Went We- = Ww! \ 


: = nt 
ce x 5" - at Ue) = Ue | 
N N tc : : 
23 [Vie -Viio] - S$ [wen wee | me 


Ve, 


E, Grid Size and Stability 


The geometry of an equilateral triangle presents a problem in 
selecting appropriate values of AY and AZ. The requirement that the 
network intersections fall on a boundary was considered essential for 
numerical calculation, dictating the relationship between AY and 4Z 
such that aY =V3«at/Z2. 

The value of AX, while arbitrary, is constrained by the fete 
that the number of X=-wise stations must be chosen such that the number 
of stations times the value of AX be equal to the predicted entrance 
length, or AX X I=X,, where I is some integer. Finally, because of 
the limited storage available in a computer, the number of grid points 
is limiced, For the purposes of this investigation, the values adopted 
were AX=.C0159, oY=0.17321, and 420.3, This grid resulted in 
25 stations in the X=direction, 11 in the Y and 6 in the Z, as shown 
in Figure l, 

Having selected the grid size, the stability of equations 13-16 
may be determined using the results of Chorin (1), who points out that 
the introduction of an artificial density results in an artificial 


speed of sound of the form: 


a 
@ ~ / § (17) 


A Mach number relative to the speed of sound is given by, 


R 2 2 2 \@ 
M= & (Unae* Vers * Wenn} (18) 
For convergence it is necessary that the Mach number (M) be less than 


one, With this condition satisfied, the system of equations is stable 


when? 
Z a x ( MIENINGUM Ne 
At ma nz (1+ 5%) DIMENSION 


(19) 
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In the present work these values were taken to be M=0.975, § =9,38025, 
4t=0,00035, and R=1.0. 
F. Friction Factor 
The wall shear stress may be related to the local pressure gradient, 
which is of importance to engineers in predicting pressure drop. The 


average stress at any section 1s given by the Newtonian shearing law, 


Fp(Se) 


Qu ; ; 
where (Sn). represents the average velocity gradient evaluated at 


(20) 


the wall PS 


Using equation 20, a dimensionless friction factor can be defined, 


7 
f= oe 





nl- 


20) 
in which f is the usual flow friction factor, 
Ge Numerical Solution 

In the present analysis, using the values of AX, AY, aZ, at 
and §S given above, equations 13-16 were found to converge to a steady 
solution after 2,000 iterations. While Chorin points out that initial 
values of variables may be set at zero throughout, the convergence can 
be hastened by using a good "guess" at the values to be calculated, 
In the present case the fully established velocity values were used 
throughout as an initial assumption, 

A copy of the program written in Fortran IV language is given in 


the Appendix. 


ig 


III, RESULTS 


The dimensionless centerplane velocity profiles are given in 
Figure 2 for five values of X/Xg (0.2, 0.4, 0.6, 0.8, and 1.0). Figures 
3 through 7 show the centerplane velocity profiles and velocity contours 
for X/X,=0.2, 0.4, 0.6, 0.8, and 1.0 respectively, The secondary ve~- 
locities, V and W, are shown in Figure 8 at a station about midway down 
the entrance region at X/X,=0.6, The value of friction factor at each 
x-wise station is plotted in Figure 9 in the form of f/fe versus X/X,, 
where the e subscript denotes the fully established values, 

The dashed line on the friction factor curve was computed from 
the local skin friction predicted from a boundary layer analysis. The 
value of skin friction must of course be unbounded at X=0 because of 
the discontinuity in velocity at the lip of the duct. For large values 
of X/Xe, f/f, must approach unity asymptotically as the results indi- 


cates 
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IV. CONCLUSIONS 


From the results of this study, it may be concluded that: 

1. The relaxation procedure used in this analysis was found to 
be convergent subject to certain stability parameters as Chorin (1) 
points out, The results, however, remain unchecked experimentally. 

2. The use of the fully established values as exit boundary 
conditions forces the relaxation values to converge on these values 
in the limit as the number of iterations increase, Therefore, the 
values of the computed velocity profiles can be only as accurate as 
the exit boundary conditions. Due to the lack of analytic and experi- 
mental results in the area, there was little choice in using the 
results of Knudsen and Katz (3) for the fully established velocity 
profile and Kays (5) for the fully established pressure gradient, 

3. In general, the grid sizes used should be reduced until the 
number of points reaches the limit of computer storage, This serves 
to maximize the accuracy of the results consistent with the accuracy 
of the fully established profiles used, 

4, It seemed intuitively obvious when first using Chorin's method 
that because it represents a convergent numerical solution to the 
finite difference form of the Navier Stokes and continuity equations, 
using an entrance length greater than a predicted value would lead to 
the fully established profile upstream of the boundary condition, 
This, however, was not the case, Thus, the necessity of using the 
predicted value of X, reduces the generality of this method in the 
solution of physical problems, The use of a Reynolds number of unity 


may have been the cause of this, however it would seem evident that a 
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complete analysis of the equations must be made numerically to find 
the reason for this limitation, 

5. The Reynolds number dependency of the governing equations can 
be extracted by using non=dimensional space parameters of the form 
X=xR/d. An analysis of the resulting equations would seem beneficial 
as the results would be completely general and would represent "uni- 


versal" profiles. 
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FIGURE 3. 


VELOCITY PROFILE AND CONTOURS 
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VELOCITY PROFILE AND CONTOURS 
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VELOCITY PROFILE AND CONTOURS 
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FIGURE 9. 
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